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We describe the tttt S-wave in D+ — » 7r+7r - 7r + decays using a unitary model for the 
■kit Final State Interactions (FSI). The three body decay is treated as a quasi two-body 
process where, at the weak vertex, the D meson decays into a resonance and a pion. 
The weak part of the decay amplitude is evaluated using the effective weak Hamiltonian 
within the factorization approximation. 

1. Introduction 

The E791 collaboration found a strong evidence for a light and broad scalar-isoscalar 
resonance in D + — > 7r + 7r~7r + decays, known as the er(500) (or /n(600)) meson 
According to their fit, this resonance contributes to approximately half of the decays 
and has a mass and width given by m CT = 478^23 ±17 MeV and = 328^40 ± 
21 MeV. This presence of the cr(500) in D + — > 7r + 7r~7r + has been confirmed by 
the CLEO collaboration in a very recent paper ^ and the results they obtained 
are in agreement with those of E791. On the other hand, this resonance was for a 
long time hidden in scattering experiments, but at present its pole position in the 
complex energy plane is rather well determined^!, which gives support to the E791 
findings. 

In our approach, we calculate the weak decay D + — > R tt + where R = <r(500), 
p(700), /o(980), using the effective weak Hamiltonian within the factorization ap- 
proximation. The (7T7r)s-wave final state is then constructed using the sc alar form 
factor previously introduced in the context of B decays to pions and kaons^El. This 
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model for the FSI, based on the framework developed in®, is unitary and takes into 
account the coupling to the KK channel. The -P-wave, which mainly corresponds to 
the p(770), and the D-wave, arising from the /2(1270), are described by the usual 
Breit-Wigner functions. The general idea of the model is depicted in Fig. 1: the 7r + 
produced at the weak vertex acts as a spectator of the rescattering between the 
other two pions. The amplitude is to be symmetrized since there are two identical 
pions in the final state. 




Fig. 1. D + decay to tt + tt ir + intermediated by a resonance R. The amplitude is symmetrized. 

2. Weak amplitudes 

To evaluate the weak decays D + — ► (<r, p, /o)vr + , we use an effective weak Hamil- 
tonian within the QCD factorization approximation. In the present case, the con- 
tributions from strong penguins are CKM suppressed and can be safely neglected 
which allows one to work at tree level. Assuming that the cr(500) is scalar-isoscalar 
with a (uu + dd) quark content, the result for the decay D + — ► (T7r + can be cast as 

As = (<nr + \H cS \D + ) = ^V cd V: d ai fAm 2 D - m*)J^(m*) , (1) 

where /„• is the pion decay constant, wid, fn-rr and m a are the D meson, pion and 
<7 masses, respectively, F®^" 7 is a transition form factor treated as a parameter 
of the model, Gf is the Fermi constant and V qiq2 is the CKM matrix element for 
the transition q\ — > (72 ■ The coefficient a± is a combination of Wilson coefficients, 
C\ + CijN c . It is worth stressing that this expression is simplified by the fact 
it receives no contribution from the color suppressed diagram since the matrix 
element (c|j M |0) with j M = 97^(1 — 7 5 )g vanishes by symmetry. The result for the 
weak decay D + — > /o(980)7r + is very similar but one must also add the non-zero ss 
component of the /o(980), which can be done introducing a mixing angle 9 in the 
spirit of ^ . The D + — > p(770) ir + amplitude does receive a contribution from the 
color suppressed diagram, proportional to 02 = C2 + C\/N c , since now (p|j p |0) 7^ 0. 

3. Final state interactions 

One can construct the 37r final state using a form factor to describe the propagation 
and decay into two pions of the resonance R. Usually a Breit-Wigner function is 
used for this purpose . In this work, we employ a model for the 5*-wave in which 



1,2008 9:19 WSPC/INSTRUCTION FILE ProceeXHadronDRB 



Resonances and weak interactions in D + — > -k + tt 7r + decays. 3 

four scalar form factors are introduced to describe the (tttt)s and (KK)s FSI SEI 
However, in the case being considered here, only one of them, namely r™*(s), en- 
ters in the calculations. This model respects unitarity and the coupling to the K K 
channel is properly accounted for. Employing an approximation where all the par- 
ticles are taken to be on-shell, this form factor is given analytically^ it is written 
in terms of the scalar-isoscalar scattering phase shifts and inelasticities^ 5 T7r (s), 
8 K x(s) and r?(s), as well as of the production functions R r {{s) and R r 2 L (s) obtained 
in chiral perturbation theory Below the KK threshold Watson's theorem is ful- 
filled: the phase of the form factor r™*(s) is given by S„„(s). The full amplitude 
D + — > <T7r + , cr — > Tt + n~ can then be written as 

Ms(u,t)=A sX [rrW + rrW] , (2) 

where As is given in Eq. ([1]). This amplitude is symmetrized as shown in Fig. 1, 
it is expressed in terms of two Lorentz invariant quantities namely u = (j>\ + 
p^-) 2 and t — (p2 + Ptt-) 2 , following the labels used in Fig. 1. In Eq. ((2]), x is a 
normalization constant related to the coupling between the pions and the scalar 
resonance, estimated to be ^1% « 30 GeV -1 . The amplitude for the decay D + — > 
/9°7r + , p° — > 7r _ 7r + is written similarly 

M P (u, t)=A P [(t- s)T p ^{u) + (u - s)T p ^{t) } , (3) 

where s = (pi +P2) 2 , r p7r7r is a relativistic Breit-Wigner^and Ap the corresponding 
weak amplitude. The D-wave is included in the model by means of another Breit- 
Wigner function and two additional parameters: a magnitude a/ 2 and a phase Sf 2 . 
The total amplitude is written as the sum of all the partial waves 

M Total (u, t) = M s (u, t) + M P (u, t) + M D (u, t) . (4) 

This expression is in a form suited to obtain a Dalitz plot in terms of u and t. 

4. Preliminary results 

Eq. (|4|) enables one to perform a fit to the E791 data reproduced as in ^21. We 
show in Tab. 1 preliminary results for the normalization constant x an d for the 
transition form factor F®~* a . The value for x is in agreement with the results from 
whereas our value for F ( f , ^°'(rnJ) is smaller than 0.79 ± 0.15 obtained in^. 



Table 1. Preliminary results for 
F^ - " 7 (m^) and for the normal- 
ization constant \- 



Par am. 


Result 




0.43 ± 0.03 


X 


22.5 ±0.8 GeV -1 
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In Tab. 2 we show our results for the fit fractions (f.f.s.), which measure the 
contribution of each channel to the total Dalitz plot. Due to interference effects, 
the sum of the f.f.s. is not necessarily 100%. In our model, the (7T7r),s fit fraction 
approximatively corresponds to the sum of the f.f.s. of a, NR and /gS of CLEO 
Compared to Oiler's work, we do not introduce a non-resonant background. The 
description of the S-wave, unified as it is in the T™*(s) form factor, avoids the use 
of a sum of Breit-Wigner functions. Use of a unitary model for the tttt FSI enables 
us to reproduce the S'-wave two-pion resonance spectrum in D + — > 7r + 7r~7r + . 



Table 2. Preliminary fit fractions from this work compared with E791, CLEO 
(isobar model fit) and J. A. Oiler results. NR stands for non resonant. 



Channel 


E79l!I(%) 


CLEcEl (%) 


Oiler 10 (%) 


This work(%) 


(7 


46.3± 9.2 


41. 8± 1.4 ± 2.5 






NR 


7.8 ± 6.6 


<3.5 


17 




/o(980) 


6.2± 1.4 


4.1± 0.9 ± 0.3 






/o(1370) 


2.3 ±1.7 


2.6 ±1.8 ±0.6 






/o(1500) 




3.4 ±1.0 ± 0.8 






(nir)s 






102 


53 


p(770) 


33.6±3.9 


20.0±2.3 ±0.9 


36 


32 


/ 2 (1270) 


19.4±2.5 


18.2±2.6 ± 0.7 


21 


10 


p(1440) 


0.7±0.8 


<2.4 


1 


2 
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